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Abstract 

We present a class of numerical solutions to the SU(2) nonlinear <r-model 
coupled to the Einstein equations with cosmological constant A > in spher- 
ical symmetry. These solutions are characterized by the presence of a regular 
static region which includes a center of symmetry. They are parameterized 
by a dimensionless "coupling constant" /?, the sign of the cosmological con- 
stant, and an integer "excitation number" n. The phenomenology we find is 
compared to the corresponding solutions found for the Einstein- Yang-Mills 
(EYM) equations with positive A (EYMA). If we choose A positive and fix n, 
we find a family of static spacetimes with a Killing horizon for < f3 < (3 ma x ■ 
As a limiting solution for (3 = (3 max we find a globally static spacetime with 
A = 0, the lowest excitation being the Einstein static universe. To interpret 
the physical significance of the Killing horizon in the cosmological context, we 
apply the concept of a trapping horizon as formulated by Hayward. For small 
values of (3 an asymptotically de Sitter dynamic region contains the static 
region within a Killing horizon of cosmological type. For strong coupling the 
static region contains an "eternal cosmological black hole". 



I. INTRODUCTION 

The aim of this paper is to discuss the static spherically symmetric solutions of the SU(2) 
nonlinear cr-model coupled to the Einstein equations with cosmological constant A > 
(subsequently referred to as the Eo"s;y( 2 )A-model). The existence of such solutions was 
suggested by the fact, that the Einstein universe || and the de Sitter spacetime || admit the 
existence of a discrete one-parameter family of static, spherically symmetric regular solutions 
of the SU (2)-0--model. These solutions are unstable, the number of unstable modes equals 
the "excitation index" . The static solutions on the de Sitter background correspond to the 
uncoupled limit, i.e. vanishing coupling constant of the Ecr,sc/(2)A-model. 

Globally regular stationary "solitonic" solutions, if they exist, play a fundamental role 
in the dynamics of gravitating systems. If they are stable, they provide possible end states 
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of evolution - and models for "stars" or "particles" . Interest in solitonic solutions and the 
subtleties of the interaction of matter with gravity was revived by the surprising numerical 
discovery of static solutions to the EYM system by Bartnik and McKinnon [H. It was soon 
found that these solutions are unstable ||, but there is a link between this instability and 
another surprise in GR: 

In recent years the discovery of critical phenomena at the threshold of black hole forma- 
tion has introduced a new twist into the gravitational collapse problem. In the original work 
of Choptuik the intermediate attractor that separates the two possible generic end-states of 
a free scalar field - the formation of a black hole and complete dispersion - is the self-similar 
"choptuon" which forms a naked singularity However, as was discovered by Choptuik, 
Chmaj and Bizon , also unstable solitons can serve as an intermediate attractor in critical 
collapse. In particular they found that the first Bartnik-McKinnon excitation actually forms 
an intermediate attractor associated with type I critical collapse phenomena - as opposed 
to the type II originally found by Choptuik (both named because of the correspondence to 
phenomena of statistical physics, for a recent overview on critical collapse phenomena see 
e.g. Choptuik @). 

While both the free scalar field and the EYM-system do not contain any dimension- 
less parameters, the a-models are simple systems with a dimensionless coupling constant. 
Thus they provide a convenient family of theories that is suitable to study genericity and 
bifurcation phenomena. 

It is known that, apart from trivial solutions, the nonlinear a- models do not admit soliton 
solutions in Minkowski space, and when coupled to gravity, there exist neither solitons nor 
static black hole solutions that are asymptotically flat (see e.g. The presence of a 

positive cosmological constant changes this situation by introducing a length scale into the 
model. From dimensional analysis one concludes that the behavior of the solutions depends 
nonperturbatively on A, and that only the sign of A is significant. 

Within the static spacetimes we restrict ourselves to spacetimes which possess a static 
region that has a regular center and is bounded by at most one Killing horizon (as opposed 
to two, such as in the Schwarzschild-de Sitter case). 

In the coupled case we find four qualitatively different types of solutions along each 
"branch" defined by the "excitation index". For small coupling constant of the cx-model, 
the solutions are globally regular, and a Killing horizon separates a static region from an 
asymptotically de Sitter dynamic region. For intermediate values of the coupling constant 
the situation is similar, but the region exterior to the cosmological horizon expands for some 
time, but eventually recollapses. For even higher values of the coupling constant the horizon 
encloses a dynamic region that undergoes collapse. Thus in some sense the static region gets 
turned "inside out". 

Solutions with a finite regular static region and positive cosmological constant only exist 
for a coupling constant smaller than some critical value which depends on the excitation 
index. The limit is singular, but combining it with a limit A — > yields solutions which are 
globally regular and static and have spatial topology S 3 . 

For A > the "first static excited states" exhibit one unstable mode, so that our solutions 
form a one parameter-family of candidate solutions for type I intermediate attractors. The 
full linear stability analysis, that leads to this result, will be presented in a separate paper 
0. 
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Similar results have been found for the EYMA system, which has been studied in great 
detail by Volkov et. al. [Q. 

The organization of this paper is as follows: In Sec. we introduce our nonlinear o- 
model and discuss its basic properties. We then specialize to spherical symmetry using the 
hedgehog ansatz for the >S£7(2)-valued matter field and write out the spherically symmetric 
field equations in a suitable gauge. 

The problem of finding regular static solutions of the Easu(2)^ problem is discussed in 
Sec. |TI[ The static regions containing the soliton are confined within a Killing horizon, out- 
side of which the timelike Killing vector becomes spacelike and the spacetimes are dynamical. 
The static regions are constructed by solving a boundary value problem as described in Sec. 
1 1 fD]. The global structure of the spacetimes is then determined by evolving the static data 
beyond the horizon. This evolution problem is discussed in Sec. [111K| . The phenomenology 
of the solutions we find is presented in Sec. [TV. 

Finally, Sec. |V| gives a discussion of our results and compares them to the results of 
Volkov et. al. for the EYMA-system. 

Conventions are chosen as follows: spacetime indices are greek letters, SU(2) indices are 
uppercase latin letters, the spacetime signature is (—,+,+,-1-), the Ricci tensor is defined 
as TZ^ V = TZ l _ l \u X , an d the speed of light is set to unity, c = 1. 



II. THE SU(2)-ff-MODEL IN SPHERICAL SYMMETRY 
A. Harmonic Maps as Matter Models 

Nonlinear a-models are special cases of harmonic maps from a spacetime (M, g^ v ) into 



some target manifold (N, Gab) (See e.g. |TT[). For the SU(2)-cr-model, the target manifold 



is taken as S 3 with Gab the "round" metric of constant curvature. Harmonic maps X A (x' 1 ) 
are extrema of the simple geometric action 



S = i[ d m x^/W\g^d,X A d u X B G AB (X). (2.1) 
Jm 



Their importance for physics was pointed out in a review article by Misner | 12| 



Variation of the matter action ( j2.1| ) with respect to the metric g^ v yields the stress- 
energy-tensor of the harmonic map 

T^u = 7 (V / ,Y a V,Y b G A b(Y c ') - l -g^X A V a X B G A B{X c )). (2.2) 

This stress-energy-tensor obeys the weak, strong and dominant energy conditions. 

In order to produce static configurations attractive and repulsive forces have to be bal- 
anced. In particular in flat space the virial theorem implies that the components of the 
stress-energy tensor cannot have a fixed sign (see the discussion by Gibbons (T]|). For 
static configurations of any harmonic map the sum of the principal pressures J2l=i ^il ls 
nonpositive everywhere |jTB[ (7« denote the spatial components of the stress-energy tensor 
with respect to an orthonormal frame). The self-interaction of any a-model can therefore 
be interpreted as attractive. In contrast for the Yang Mills field the sum of the principal 
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pressures is nonnegative, which can be interpreted as a repulsive self-interaction. Both fields 
therefore do not allow soliton solutions on a flat background. But while the Yang-Mills field 



coupled to gravity admits solitons by cancelling the repulsive force with gravity [0, the 
gravitating a-model correspondingly does not (see e.g. [(J). 

Thus we cannot expect static solutions to exist unless we add some "repulsive force" - 
or in the presence of nontrivial topology. As an example for the latter one could think of 
a static spherical universe of topology S 3 , where there is a balance between the tendencies 
to collapse towards either "center" . A repulsive force on the other hand can be introduced 
with a positive cosmological constant A. Therefore we consider the total action 

S = [ d A x^{-^-{Tl - 2A) + C M ), (2.3) 
/ 1o7tG 



where Cm is given by the Lagrangian of (|2.1|) . This action gives rise to the Einstein equations 

- -g^K + kg^ = 8nGT^, (2.4) 

as well as to the field equations of the harmonic map 

<T (V M V,X A + Ti c (X D )V,X B V u X c ) = 0, (2.5) 

where T^ c denote the Christoffel symbols with respect to Gab- The field equations thus 
allow a general nonlinear dependence on the fields X A through the Christoffel symbols and 
they depend quadratically on the first derivatives of the fields. 

In units where c = 1 the coupling constant 7 of the harmonic map has dimension 
mass /length, whereas the gravitational constant G is of dimension length /mass. Both 
constants enter the equations only in the dimensionless product (3 = AnG'-f, thereby defin- 
ing a one-parameter family of distinct gravitating matter models. The parameter A plays 
a different role, since it has dimension 1/ length 2 . Thus, when the cosmological constant 
is nonzero, it provides the length scale of these theories. Therefore all theories with the 
same value of (3 are equivalent irrespective of the value of A. If, on the other hand, A = 0, 
the field equations are scale invariant. The gravitating SU(2)-<r-model with non-negative 
cosmological constant thus corresponds to a two-parameter family of inequivalent theories, 
parameterized by the continuous parameter (3 and the discrete parameter sign(A). 



B. Spherical Symmetry 

A spherically symmetric metric can be written in the general form as the warped product 
(see e.g. JT§) 

ds 2 = dr 2 + R 2 {p,t)dQ 2 , (2.6) 

where dr 2 is a general two-dimensional Lorentzian line-element and the area of the orbits 
of SO(3) is given by 4nR 2 (p,t). Choosing the coordinates t and p orthogonal one can write 
the metric in the form 
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ds 2 = -A(p, t)dt 2 + B(p, t)dp 2 + R 2 (p, t)dtt 2 . 



(2.7) 



In ( [2.7|) there is some gauge freedom left, which can be used to eliminate one of the three 
functions A, B or R. 

One way to fix the gauge would be to use the function R as a coordinate. This is possible 
as long as V M i? 7^ 0. This gauge is usually referred to as the Schwarzschild gauge. Here we 
will deal with the situation that V^i? becomes zero on some maximal two-sphere (which is 
not a horizon), a phenomenon which we will discuss in detail below. 

Following Volkov et al. Jj]] we therefore choose a different gauge, which keeps close to the 
standard Schwarzschild line element in another way: we fix A(p,t) = 1/B(p,t) = Q(p,t) 
and keep the area of the SO (3) orbits as the second free function: 

ds 2 = -Q(p, t)dt 2 + — dp 2 + R 2 (p, t)dQ 2 . (2.8) 

Q(p,t) 

The metric ( |2.8| ) is regular where Q(p,t) and R(p,t) are regular functions on spacetime, 
except possibly at points where either Q(p,t) or R(p,t) vanish. The vanishing of Q(p,t) is 
related to the existence of horizons, which will also be discussed below. 

For a regular spacetime the vanishing of R(p, t) is associated with the singularity of 
spherical coordinates at the axis of symmetry. Such a coordinate center does not necessarily 
have to be present (e.g. in a spherical wormhole), but we will require our spacetimes to 
possess at least one regular center - there may also be two, such as in the Einstein static 
universe written in spherical coordinates. In the following, several boundary conditions will 
be derived from regularity requirements, we will therefore make some comments on regularity 
near the center of spherical symmetry for the metric defined in Eq. (|2.8j ) : We assume the 
existence of four regular (meaning C°° or C k as appropriate) coordinate functions x, y, z, 
t on the manifold. All other functions are defined as regular if they can be expressed as 
regular functions of x, y, z, t. Furthermore, the relations between the functions R, 8, 
and x, y, z are required to be the standard coordinate transformation between Cartesian 
and spherical coordinates: x = .Rsin^sin^, y = R cos (j) sin 8 , z = Rcos8. Note that the 
spherical coordinates, in particular R, are not regular functions, but all even powers of R 
are. Near the axis we choose the parameterization of the radial coordinate p such that it 
has the same regularity features as R: 

p = Rh(R 2 ), 

where h is a regular positive function - and thus also p itself is not regular. Any spherically 
symmetric function which is a regular function of x, y, z, t (or just "regular" for short) can 
therefore be written as a regular function of p 2 and t. 

For any fixed t the center R = then is a regular point of spacetime, if Q(0, t) = h(0). 
Without restricting generality we can choose h(0) = 1 for convenience, such that 

R(p,t)=p(l + 0(p 2 )) 

Q(p,t) = l + 0(p 2 ). (2.9) 
In spherical symmetry it is possible to define a quasilocal mass 
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m = -{l-V li KW l R), 



(2.10) 



which reads 

m = f( 1 -f -Q(R') 2 )- (2-n) 

in our coordinate system (|2.8|) . For a recent discussion of the properties of the quasilocal 
mass in spherical symmetry see Hayward [ |T6|j . 

The harmonic map field configuration X A (x^) can be called spherically symmetric, if 
the Lie derivatives of the energy-momentum tensor with respect to the Killing vector fields 
that generate the SO (3) action vanish. One possibility to achieve this, is to demand that all 
fields be functions of p and t only. This would leave us with a coupled system of differential 
equations for three fields and two metric functions. Since the target manifold (S 3 , G) also 
admits SO (3) as an isometry group, there is an alternative way to impose spherical symmetry 
on the harmonic map, which is the well known hedgehog ansatz. We first introduce spherical 
coordinates (/, O, $) on the target manifold, writing the SU(2) line element as 

ds 2 = df + sin 2 f{dG 2 + sin 2 6rf$ 2 ). (2.12) 

The hedgehog ansatz now ties the coordinates on the target manifold to those on the base 
manifold: 



f(x») = f(p,t), Q(x») = 0, (2.13) 

Due to this ansatz two of the three coupled fields are already determined and only one field 
f(p,t) enters the equations. The matter field equations ( |2.5p are then reduced to the single 
nonlinear wave equation 

□/=^, (2.14) 

where □ is the wave operator of the spacetime metric. 

Regularity of the harmonic map within this ansatz means that when expressing the field 
in terms of regular coordinates X A on SU(2), the X A are regular fields on spacetime. Since 
X {x^) = (f(p,t)/p)x A , where x A denote the Cartesian coordinates x,y,z on spacetime, 
the field f(p, t) thus has to be of the form 

f(p,t)=pH(p 2 ) = f'(0)p + O(p 3 ) } (2.15) 

where if is a regular function. Thus / vanishes at the center, which means that the center 
R = is mapped to one of the poles of the target manifold, that have been fixed via the 
hedgehog ansatz. 
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III. STATIC SOLUTIONS 



A. Static Field Equations 

In addition to spherical symmetry we assume that spacetime admits a hypersurface 
orthogonal Killing vector field d t , which is timelike in some neighborhood of the (regular) 
center R = 0. 

The combinations ( f t ) — (£) — 2(g) and (\) — (?) of the mixed components of Einstein's 
equations plus the matter field equation ( |2.14| ) then yield the following set of coupled second 



order autonomous ODEs for the metric functions Q(p),R(p) and the matter field f(p): 

(R 2 Q')' = -2AR 2 , (3.1) 
R" = -(3Rf 2 , (3.2) 
(QR 2 f')' = sin 2f. (3.3) 

Furthermore the above system of equations has a constant of motion: 

2(3 sin 2 / + R 2 (A - (3Qf 2 ) + RQ'R' + QR' 2 -1 = 0. (3.4) 

This expression can either be derived by integrating the system (fj . (PT5|) and using the 
regularity conditions at the axis, or by using the ( p p ) component of Einstein's equations. 

Note that if / is a matter field solution, then so is also k 7r ± / for any integer k. 

If A is nonzero it sets the length scale, it thus can be eliminated in these equations 
by using the dimensionless quantities p = VXp and R = \/~K.R. If, on the other hand, A 
vanishes, the equations are scale invariant, and thus invariant under rescalings p = ap and 
R = aR. Furthermore, in this case Eq. (|3.1|) together with regularity conditions at the axis 
implies Q = 1. 

If we integrate the equations with the above boundary conditions from p = to larger 
values of p one of the following four situations has to occur: 

1. the static region "ends" in a singularity, 

2. integration might run into a second (regular) pole R = 0, which would mean that the 
resulting spacetime has compact slices of constant t and is globally static, 

3. the static region might persist up to spatial infinity, 

4. the static region might be surrounded by a Killing horizon beyond which spacetime 
becomes dynamical. 

The first case can easily be produced by shooting off a regular center with arbitrary initial 
data, it can be excluded by setting up a boundary value problem that enforces one of the 
other three cases. The second case can be discarded for our model with positive cosmological 
constant, since the existence of a static region with two regular centers is not compatible 
with the field equations: To see this recast Eq. ( |3.1| ) into the integral form ( |3.21| ). It is clear 
then, that - for nonzero A - Q' diverges if R goes to zero a second time. For A = of course, 
such solutions may exist as e.g. the static Einstein universe (|3.7|), which will be discussed 
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in the next section. The third case can also be excluded if the cosmological constant is 
positive: a static region can not be extended to spatial infinity, but rather a singular point 
of the equations Q(ph) — 0, that corresponds to a horizon, has to develop. This can be seen 
as follows: Suppose the solution exists up to p — > oo and R is monotonically increasing - all 
other assumptions automatically lead to one of the cases 1, 2 or 4 - then one can show, that 
for A positive Q' < — const./ p for p large enough, which means that Q would be bounded 
from above by a function that tends to — oo as p —>■ oo. So again Q has to cross zero at 
some finite value of p. 

For positive cosmological constant we therefore may confine ourselves to the cases with 
horizon, and we will construct static regions as an ODE boundary value problem, where 
the boundary conditions correspond to a regular center at p = and a regular horizon at 
p = ph, where the determination of the value of pn is part of the boundary value problem. 
The appropriate boundary conditions at the horizon will be determined in Sec. [Ill C[ and 
the boundary value problem is described in Sec. [I I HI]. 

Given a static region, i.e. a region of spacetime where the Killing vector d t is timelike, 
which is bounded by a Killing horizon, we can consider as a second step the time evolution 
problem of the static data on the horizon into the dynamic region where the Killing vector 
d t is spacelike and the spacetime is thus homogeneous. The time evolution problem thus 
reduces to a system of ODEs which are solved as an initial value problem as described in 



Sec. HIE. 



B. Exact Solutions 

Some solutions of Eqs. with a regular center can be given in closed form. 

They do arise as certain limits of the numerically constructed family of solutions to be given 
below. 

First of all, for A > 0, for the trivial case / = Eq. (|3.2|) has the solutions R(p) = ap + b. 
Imposing the regularity conditions ( |2.9| ) this gives de Sitter spacetime 

R(p)=p, Q(p) = l-±£, f = 0. (3.5) 
For / = 7r/2, R(p) is of the same form as above. For b = we get 

R(p)=p, Q(p) = 1-2(3 -±f, (3.6) 

which looks like de Sitter for large p but has a conical singularity at the center if (3 > 0. 
Furthermore the static region shrinks with increasing (3 and ceases to exist for (3 = 1/2. In 
the limit of vanishing coupling constant (3 = 0, where spacetime is de Sitter, this solution 
/ = 7r/2 is still singular at the center with diverging energy density. Nevertheless the total 
energy is finite, and this solution may be viewed as the "high excitation" limit of the regular 
solutions, that exist on de Sitter background ||. 

Finally, for A = Eq. ( |3.1| ), together with regularity conditions at the axis, yields 
Q = 1. For (3 = 1 the remaining equations can be solved analytically to give the static 
Einstein universe: 
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R(p) = ship, Q(p) = 1, f(p) = p. 



(3.7) 



Note that the stress-energy tensor, which has the form of a perfect fluid in this case satisfies 
p + 3p = 0. As will be described in detail in Sec. [VC the static Einstein universe arises 
in the limit of maximal coupling constant of the numerical constructed first excitation iff at 
the same time A is set to zero. 

For completeness we mention that there also exist exact solutions to Eqs. (p.lp-flOp. 
that do not posses a regular center of spherical symmetry, such as the Nariai spacetime. 



C. Horizons and global structure 

In order to discuss the global structure of the spacetime and in particular regularity 
questions from which we derive the boundary conditions for our ODEs, it is helpful to 
consider a coordinate system, which is regular at the horizon. We write the metric (|2.8|) as 



ds 2 = -Q{p)du 2 - 2du dp + R(p) 2 dtt 2 , (3.8) 

where the coordinate function p and the metric function Q(p) coincide with those in the 
metric (|2.8|) , and the coordinate u is given as 



Note that the coordinates ( |3.8| ) cover only half of the maximally extended spacetime. 
In the following, we will simplify our discussion by only talking about the Killing horizon 
contained in the portion of spacetime covered here. All statements made can be extended 
trivially to the complete spacetime and in particular the second component of the horizon 
by time reflection. We also remark that all solutions have the topology S 3 x R. 

The static Killing vector field is d/du = d/dt, where the latter is taken with respect to 
the (t, p) coordinates. The metric fl3.8|) is regular if Q(p) and R(p) are regular functions, 
except when R = 0, which corresponds either to the usual coordinate singularity of spherical 



symmetry, which has been discussed in Sec. (|II B| ) or to a spacetime singularity, as discussed 
in Sec. ( PIE|) . 

The Killing vector field d u , which we have assumed to be timelike in some neighborhood 
of the (regular) center R = need not necessarily be globally timelike. Regions where it 
becomes spacelike, i.e. Q(p) < 0, are dynamic with homogeneous spacelike slices of constant 
time p, the spatial topology being S 2 x R. Such regions thus correspond to Kantowski-Sachs 
models. At the boundary of static and dynamic regions the metric function Q(p) vanishes, 
such a surface of p = const, is thus a null surface. Furthermore the Killing vector field d u is 
null and tangent on this surface, which therefore is a Killing horizon. 

In the presence of an asymptotic infinity, such as in an asymptotically flat or asymptot- 
ically de Sitter spacetime, one can use the asymptotic region to classify event horizons, e.g. 
as black hole or cosmological event horizons. Furthermore one is provided a straightforward 
definition of "inward" and "outward" directions. In the cosmological case these issues are 
less clear. We therefore follow the definitions of Hayward [Qj: He proposed a local defini- 
tion of a trapping horizon, a concept which does not make use of asymptotic flatness and 
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is therefore also suitable for more general situations. Intuitively, the physical interpretation 
of the Killing horizon depends on whether the dynamical region is collapsing or expanding 
off the horizon, and whether this region is to be interpreted as inside or outside. Both can 
be formulated in terms of the null expansions Q± of ingoing and outgoing null rays from 
two-surfaces. In spherical symmetry we consider R = const, surfaces, with null expansions 

e± = ^£±R\ (3-io) 

where C± is the Lie-derivative along the null directions 

l+ = d p and l^ = 2d u -Qd p (3.11) 



respectively, so 



0+ = 2% and 0_ = -2Q^-. (3.12) 
R R 



as a 



The first use of the expansions is to define a trapped surface in the sense of Penrose |L7 
compact spatial surface for which B_0 + > 0. If one of the expansions vanishes, the surface 
is called a marginal surface. For a non-trapped surface 0_ and + have opposite signs, and 
we call directions in which the expansion is positive outward, and inward when it is negative. 

On the Killing horizon Q = 0, 0_|q =o = while + |q=o ^ and £ + 0_|q =o ^ (except 
when also R' = 0, which we exclude for the moment). In Hayward's terminology || such 
a three-surface is called a trapping horizon. It is said to be outer if £ + 0_|q =o < 0, inner 
if £ + 0_|q =o < 0, future if + |q =o < and past if + |q =o > 0. A future outer trapping 
horizon provides a general definition of a black hole, while a white hole has a past outer 
horizon, and cosmological horizons are inner horizons. 

If R'(p) = within the static region then 

£+0-|q=o<O, 0+|q=o<O (3.13) 

and the surface Q = is a future outer trapping horizon. On the other hand if R' vanishes 
in the dynamical region or is monotonic, then the signs in ( |3.13|) are reversed and one speaks 
of a past inner trapping horizon. 

Note, that at points where R' — 0, both ± = 0. However, this does not lead to 
a (trapping) horizon since the expansions merely change sign. The significance of such 
marginal surfaces is that the meaning of inward and outward directions are reversed. A 
typical example is an "equatorial" two-surface of a round three-sphere at a moment of time 
symmetry. 

Using the expansions + and 0_ (3.12) one can rewrite the quasilocal mass ( |2.1l ) 



m = f (1 + T e+e - ) - (3 - 14) 

Using the above assignment of inward and outward direction the quasilocal mass can be 
interpreted as the total mass that is contained within any spatial three volume, that is 
bounded by the two-sphere with areal radius R. On the marginal surfaces, where R' = 0, 
the quasilocal mass equals R/2. 
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D. The boundary value problem for the static region 



The system of ODEs (|3.1| ) -(|3.3| ) has singular points for R = and Q = 0. Making use 
of the imposed regularity conditions at the axis ( |0| ) and ( |2.15| ) the equations demand in 
addition Q"(0) = — 2A/3. Note that solutions, that are regular in a neighborhood of p = 
are determined by the value of the single parameter c = f'{0). 

The other singular point occurs at the horizon, where Q{ph) = 0. Formal Taylor series 
expansions around p = pu give 

Q{p) = Q'h(p - Ph) + ^^(p - Ph? + 0((p - p H f) 

R(p) = R H + R'(p H )(p - PH ) + ^f^-{p - p H f + 0((p - p H ) 3 ) 

f(p) = b + f'(p H )(p - Ph) + A^(P " Ph) 2 + 0((p - p H f), (3.15) 

where 

p H , f(p H ) = b, R( Ph ) = Rh, Q'(ph) = Q' h (3-16) 

are free shooting parameters. The other coefficients are determined by the requirement that 
Q and R be regular functions of p. Consistency with Eqs. ( |3~T|) - (|3.4j ) then amounts to the 
conditions: 

1 - 2(3 sin 2 b-R\ A 
R [Ph) = ^-^r, , (3.17) 

Q"(p H ) = _ 2 R ' ipH l Q ' ipH) - 2A, (3.18) 



R H 

sin 2b 

Q'(ph)R 2 „ 



f'(P») = T^T • (319) 



„»/ n sin 2b , R'(ph) A . 

f (Ph) = nH i p2 ^ + cos 26 + -— - . 3.20 

Q{Ph)R h r h Q\Ph) 

In order to determine the spectrum of static solutions for a fixed coupling constant j3, we 
solve the boundary value problem Eqs. ( |3~T|) - (|3.3j ) with boundary conditions ( |2.9|) , (|2.15|) 



and ( |3.15| ) between axis and horizon. We use a standard two point shooting and matching 
method (routine d02agf of the NAG library [0), where the parameters f'(0),pH, /(.Ph), Rh 
and Q' H serve as shooting parameters. 

For (3 = a discrete one-parameter family of solutions has already been discussed in 
@]. In order to get good initial guesses for the shooting parameters for (3 > 0, we follow 
one solution from (3 = up to higher values of j3, interpolating the values of the shooting 
parameters at the present and last "/3-step" to obtain values for the next "/9-step" . 



E. Integration Through the Horizon 

Since the singularity at the horizon is merely a coordinate singularity, we can extend 
spacetime through the horizon and reintroduce the coordinates in (|2.8[) in the dynamic 
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region beyond the horizon. The Killing vector field dt becomes spacelike, and instead of t 
the timelike coordinate is p. The coupled system of ODEs ( |3.1|) - ( |3.3|) together with initial 
conditions ( |3.15| ) therefore constitute an initial value problem. 

While integrating forward in time p, essentially two things can happen according to the 
behavior of R(p): 

1. R{p) is monotonically increasing for all p > 0. Then time evolution exists for all p > 0. 
This can be seen by turning Eqs. (|3~l| )- (|3.3| ) into integral equations: 



p 

2A ' 



Q' = -— I R 2 dp (3.21) 



o 



R'=l-P J Rf' 2 dp (3.22) 



/' = ^ I sin(2/)dp. (3.23) 



o 

p 







All first derivatives are bounded as long as Q and R don't go to zero. Beyond the 
horizon Q cannot go to zero, since Q' < for all p > and Q(pu) = 0. R cannot go 
to zero since it is monotonically increasing by assumption and -R(O) = 0. Therefore 
Q, R and / stay finite for all finite p. 

From Eq. ( ggg ) it follows that < R' < 1 for all p > 0, so R = O(p) for p -»• oo. 
From Eq. Q3.21| ) we get, that Q = 0(p 2 ) and from Eq. ( p.23j ) we see that /' goes to 
zero as p~ 3 and therefore / goes to a constant at infinity. 

2. R develops an extremum at some finite p extr > 0. Since R" < for all p and R'{p) < 
for all p > p ex tr, which follows from Eq. ( |3.22| ), R goes to zero at some finite ps > pn- 



Here we already excluded the case ps < Ph in Sec. |III A| , since this means a singularity 
in the static region. 

For the same reasons as above, time evolution exists for all p < p$- In the limit p — > 
Ps a spacetime singularity occurs. This follows from inspection of the Kretschmann 
invariant: 

R^R^ar = ( 4 + 4R 2 Q ,2 R' 2 + R A Q" 2 + 

+ 8QR' (R 2 Q'R" - R') + AQ 2 (R' 4 + 2R 2 R" 2 )). (3.24) 



Since by assumption Q(p), Q'{p) and R'{p) are negative near ps all terms of Eq. (|3.24|) 



are non-negative and at least some of them clearly have a nonzero numerator while 
the denominator vanishes with some power of R. 

The construction of solutions beyond the horizon consists of solving the initial value 
problem, i.e. we integrate Eqs. (|3.1|)-(|3~3]) with initial conditions ( |3.15| ) for p > p H , where 
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the parameters Ph, f(pii), Rh and Q' H are determined by the solutions of the boundary 
value problem described in Sec. 1 1 1 D . For numerical integration we used routine d02cbf of 
the NAG library |8|. 



IV. PHENOMENOLOGY OF SOLUTIONS 

A. Phenomenology of Numerically Constructed Solutions with A > 

For (3 = the Eer^^A equations decouple into Einstein's vacuum equations with A and 
the matter field equation (|2.14j ) on the fixed background. With our regularity conditions at 



the axis the solution to Einstein's equations is de Sitter space 



= R(p)=p, (4.1) 



whereas the field equation ( |3.3|) admits a discrete one parameter family of regular solutions 
0]. Within the static region these solutions oscillate around 7r/2 where energy increases 
with the number of the oscillations and is limited from above by the energy of the "singular" 
solution / = 7r/2. Outside the cosmological horizon at p = 1 the solutions remain finite and 
tend to a constant near infinity. 

Increasing the coupling constant (3, the numerical analysis shows, that solutions of this 
type persist as long as (3 does not get too large. The qualitative behavior of the field / in 
the static region is the same as in the uncoupled case, i.e. the n-th excitation oscillates n 
times around n/2, whereas the behavior of the field in the dynamic region as well as the 
behavior of the geometry depend strongly on the value of the coupling constant (3. 

To summarize, we get the following qualitative picture of solutions in dependence on the 
coupling constant (3: 

• For small (3, < (3 < [3 crit {n), the solutions are similar to those of the uncoupled case 
(3 = 0, that is: the area of SO (3) orbits is monotonically increasing with p, beyond 
the horizon the solutions persist up to an infinite value of the coordinate time p. Near 
infinity the geometry asymptotes to the de Sitter geometry, that is R = O(p) and Q 
tends to — oo as 0(p 2 ). According to Hayward's definitions ]2[ the horizon is an inner 
past trapping horizon, separating the static region from an expanding dynamic region. 
The field / shows the same qualitative behavior as in the uncoupled case (3 = (see 
Figs. [I] and |). 

• For (3 = f3 cr it{n) the areal radius R(p) still increases for all p > but this time goes 
to a constant at infinity, that is R'(oo) = 0. For even stronger coupling, (3 cr i t {n) < 
(3 < (3*(n), R develops a maximum in the dynamic region, i.e. R'(pe) = at some 
finite time Pe > Ph (see Fig. |). As was discussed in Sec. [Ill E| R then decreases 
and goes to zero at some finite coordinate time ps ~ which corresponds to the finite 
proper time t$ = f p ^ dp/ a/ Q(p)- This causes the geometry to be singular at ps- The 
horizon, which again is an inner past trapping horizon separates the static region from 
an initially expanding dynamic region, which reaches its maximal spatial extension 
aX, pe and then recollapses to a singularity at p$. The maximum of the areal radius 
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occurs at earlier and earlier times as the coupling constant is increased until it merges 
with the location of the horizon when j3 = 



• At (3 = /3*(n) the "maximal two sphere" coincides with the horizon, pe = Ph- The 
fact that R'(ph) = at this value of (3 may be interpreted as exchanging the inward 
and outward direction at the horizon: for smaller values of f3 the static region was 
surrounded by the inner past trapping horizon, whereas for larger values of the cou- 
pling constant, < (3 < (3 maX) where pe < Ph, the static region encloses the 
horizon, which becomes now an outer future trapping horizon. Beyond this horizon, 
the dynamic region undergoes complete collapse at p = p$. 

Embedding diagrams of the static regions of the first excitation for several values of 
(3 can be found in Fig. [|. We note, that the "critical" values of the coupling constant, 
f3crit{n), (3*{n) and (3 max {n) decrease with the excitation number n. In Sec. [IV C] we will give 
an argument, that f3 max (n) is a decreasing sequence, which is bounded from below by the 
maximal coupling constant (3 max (oc) = 1/2 for the "singular" solution fl3.6|). 

The solutions described above exist in the presence of a positive cosmological constant 
A > 0. As will be described in detail in the next section, Sec. |1V B| the limit (3 — > (3 max yields 



regular solutions iff one takes the limit A — > appropriately as described in Sec. |IVB| . 



B. The Limit f3 ->■ (3 max (n) 
Recall from Sec. [Ill A| that the cosmological constant A sets the length scale in Eqs. (ETT 



and that it can be eliminated from these equations, by introducing the dimensionless 
quantities p = \fhp and R = \fhR. This corresponds to measuring all quantities that have 
dimension of length, as e.g. the energy E, the coordinate distance of the horizon pn from 
the origin, the radial geometrical distance of the horizon dn from the origin, the areal radius 
Rh of the horizon, and l//'(0), in units of l/v^A. We find that all parameters, that have 
dimension of length go to zero in the limit (3 —>■ (3 max when measured with respect to this 
length scale. This indicates that 1/y/X is not the appropriate length scale for taking this 
limit. We therefore switch to the alternative viewpoint of pn as our length scale, and we 
fix pn = 1- In this setup A depends on (3 and the excitation index n and goes to zero in 
the limit (3 — > (3 max . The parameters E, du and l//'(0) attain finite values when measured 
in units of pn, whereas Rh/ph goes to zero. (See Fig. ||). This strongly suggests, that 
there exists a solution with (3 = (3 max which obeys Eqs. (|Q|)- (|3l4T ) with A = and has two 
centers of symmetry. In particular this means that the static region of this solution has no 
boundary, since any t = const slice has topology S 3 . 

Furthermore, as can be seen from Fig. |6|, the dimensionless parameter f(pii) for the first 
excitation tends to tt, and R'(pn) tends to —1 in the limit j3 —>■ (3 max . As will be shown in 
the next section, Sec. |IV Q A = implies Q = 1. The limiting solution with A = will 
therefore satisfy the regularity conditions ( |2.9| ) and ( [2.15| ) not only at the axis p = but 
also at the second zero of R, which means that such a solution is globally regular with two 
(regular) centers of spherical symmetry. In fact, for the first excitation this limiting solution 
is just the static Einstein universe (|3.7|), which can be given in closed form. 
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These observations allow one to determine the maximal value of the coupling constant 
/3max{n) not as a limiting procedure (3 — > /3 max , but rather by solving the boundary value 
problem Eqs. ( |3.1| ) - ( |3.3| ) with A = and with boundary conditions, that correspond to 
two regular centers of symmetry. 



C. Globally Static, regular Solutions for A = 

For A = Eq. ( |3.1| ) can be solved immediately to give R 2 Q' = const. According to the 
regularity conditions at the axis (|2.9| ) the constant has to vanish, which means that Q' = 
and therefore Q = 1. The remaining system of equations is: 

R" = -(3Rf 2 , (4.2) 

(R 2 f )' = sin(2/) (4.3) 

and 

2(3 sin 2 / - (3R 2 f 2 + R' 2 - 1 = 0. (4.4) 

Note, that this system of ODEs is scale invariant, that is any solution R(p),f(p) leads via 
rescaling to the one parameter family of solutions given by aR(ap), f(ap). Keeping this in 
mind, we can fix the scale arbitrarily, e.g. in setting the first derivative of the field / equal 
to one at the origin: f'(p = 0) = 1. Thereby any solution, that is regular at the origin, is 
determined entirely by the value of the coupling constant (3. 

Regularity conditions at the second "pole" R(pp) = are the same as at the origin, except 
that / either tends to n, if its excitation number is odd, or to if it has even excitation 
number. This can be inferred from 7r/2 < f(pn) < 7r for n odd and < f(pn) < 7! ~/2 for n 
even for all (3 < (3 max , since this is the case for (3 = and according to (|3.15|) no crossing 
of the zero-line or 7r-line is allowed. Note that this corresponds to all odd solutions having 
winding number 1, whereas even solutions are in the topologically trivial sector. 

These regularity conditions together with the invariance of the equations under reflection 
at the location of the maximal two-sphere R'(pe) = 0, causes globally regular solutions R(p) 
to be symmetric around pe whereas f(p)—ir/2 is either antisymmetric for n odd or symmetric 
for n even. 

For / symmetric the formal power series expansions of R(p) and f(p) around p = Pe 
gives 

R( P ) = R(Pe) + 0((p-pe) 4 ), 

f(p) = arcsin + ^~^ ^f^ + 0((p - p E )% (4.5) 

and for / — tt/2 antisymmetric we get 

R(p) = ^jj^ 2 - (2/3 - l) (P ~r )2 + 0(( P - peY), 

f(p) = I + f(PE)(p ~ Pe) + 0((p - pe) 3 ). (4.6) 
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In order to solve the system ( f4.2|) , ( |4.3|) we again use the shooting and matching method 
on the interval [origin, pe] using the above Taylor series expansions to determine the bound- 
ary conditions at p = pe- Shooting parameters are now p E , /'(pe) and (3 for odd solutions 
and pe, R{.Pe) and (3 for even solutions. The results are displayed in Table |. 

It is clear from (|4.5|) and ( |4.6| ), that regular solutions for A = can only exist if (3 > 1/2. 
Assuming now, that our numerical observations concerning the first few excitations extend 
to higher excitations, we give the following argument: Since every "branch" of the "A > 
solutions" persists up to a maximal value of beta, which can be computed by solving the 
boundary value problem ( |4.2[ ) together with regularity conditions at the two "poles" - which 
implies (3 > 1/2 - and since we know, that in the limit (3 — > there exists an infinite number 
of excitations f§], we conclude that this whole family of solutions with A > persists up to 
some maximal value /3 ma x{n), which is greater than 1/2. In other words, for any (3 < 1/2 
there exists a countably infinite family of solutions with A > 0, whereas for (3 > 1/2 our 
numerical analysis shows, that only a finite number of solutions exists. (See Table |). 

V. DISCUSSION AND OUTLOOK 

We have shown numerically that the SU(2)-cr-model coupled to gravity with a positive 
cosmological constant admits a discrete one-parameter family of static spherically symmetric 
regular solutions. These solitonic solutions are characterized by an integer excitation number 
n. A given excitation will only exist up to a critical value of the coupling constant (3; the 
higher n, the lower the corresponding critical value. Our calculations indicate that the 
infinite tower of solitons present on a de Sitter background persists at least up to a value 
of /3 — 1/2. Thus there exists a /3 > 1/2 beyond which the number of excitations is finite 
and decreases with the strength of the coupling. As mentioned, qualitatively the cx-model 
under consideration shows striking similarities to the EYM system as studied in detail by 
Volkov et.al. The main difference being that the static solutions to the EYM-system depend 
on the value of the cosmological constant while in our case A scales out from the equations 
and (3 plays the role of a "bifurcation" parameter. Another difference concerns the globally 
regular static solutions with compact spatial slices. For the EYM system these appear for 
definite values of A(n) while for the cx-model the corresponding solutions exist only in the 
(singular) limit as A goes to zero and definite values of (3. Thus in our case there are closed 
static universes with vanishing cosmological constant, the lowest excitation being the static 
Einstein cosmos. This is possible because in this case the stress-energy tensor of the cx-field is 
of the form of a perfect fluid with the equation of state p = Another interesting aspect 

is the geometry of a given excitation as a function of the coupling strength: the static region 
is always surrounded by a Killing horizon separating the static from a dynamical region, 
which for small couplings becomes asymptotically de Sitter. As the coupling is increased 
the two-spheres of symmetry beyond the horizon are first past and then become future 
trapped and a cosmological singularity develops. Finally, for even stronger couplings, again 
the region beyond the horizon collapses, but within the static region the in- and outgoing 
directions (as defined by the sign of the expansion for null geodesies) interchange. 

An important question to be answered is whether these solitons are stable under small 
radially symmetric time dependent perturbations. In a forthcoming publication we intend 
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to present a detailed stability analysis. We will show that for A > all excitations are 
unstable with their number of unstable modes increasing with n. This was to be expected 
at least for small coupling. The lowest excitation thus has a single unstable mode and it is 
known, from other models, that such a solution can play the role of a critical solution in a 
full dynamical treatment of spherically symmetric collapse. 
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TABLES 



n 


ftmax 


PP = dp 


E/4-irj 


E/Air^dp 


1 


1 


7T 


3vr/2 


3/2 


2 


0.74255 


6.74225 


11.78039 


1.74724 


3 


0.64931 


12.10140 


22.43662 


1.85405 



TABLE I. Results for the first three excitations for A = 0. Since Q = 1 the coordinate distance 
pP of the two regular "poles" equals the radial geometrical distance dp. The energy density pp 
and energy E are given in units where /'(0) = 1. The ratio E/dp can be compared to the results 
for solutions with A > and represents the limit (3 — > f3 max for those solutions. 
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FIGURES 



FIG. 1. The first three solutions for coupling constant j3 = 0.3 < Peru- Inside the horizon at 
p/ph = 1 the "n-th" excitation crosses the 7r/2-line n times. In the dynamic region outside the 
horizon the solutions evolve towards a constant. 



FIG. 2. The energy density u n = Qn{f' n ) 2 + sin 2 (/ n )/i? 2 for the first three excitations at 
[3 = 0.3. - denotes the frist, x the second and + the third excitation, all p, n have been nor- 
malized to unity at the axis. The corresponding total energies, measured in units of 1/y/A. are 
Ex/Airy = 1.66621, E 2 /^ = 1.88082 and E 3 /4irj = 1.956613. 



FIG. 3. The area function R(p) for the first excitation for (3 near /? C rit(l)> 
0.470366 < I3 crit {l) < 0.470373. The vertical line marks the horizon at p = 0.88761. 



FIG. 4. The rotational surfaces z(R) correspond to the embedding of a t = const-slice 
(0 = tt/2) into R 3 . The upper and lower half part of the diagram resemble the two static re- 
gions causally separated by the horizon at z = 0. Shown is the geometry of the first excitation for 
(3 in the range < (3 < 0.913 < (3 max . The sphere of unit radius for de Sitter space ({3 = 0) gets 
more and more deformed as (3 increases. 



FIG. 5. Some parameters with dimension length measured in units of pn- Except Rh all of 
them stay finite in the limit (3 — > (3 max . Moreover in the limit the values tend to the corresponding 
ones of the static Einstein universe, as given in table Q. 



FIG. 6. The (dimensionless) parameters f{pn) = b and R'(ph) of the first excitation in depen- 
dence of the coupling constant (3. /(ph) tends to ir in the limit [3 — > (3 max and R'(ph) tends to — 1 
in this limit. 



20 



3 




\og(p/p H ) 
Figure 1 



Mo) 




O 0.2 0.4 0.6 0.8 1 



p/rhoH 
Figure 2 



2 



1.5 - 




0.0001 0.01 1 1e+02 1e+04 

P 



Figure 3 



-1-5 



1.5 



0.5 




R 



0.5 



1.5 



Figure 4 









/'((»/'// 













b 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

P 

Figure 5 



